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Abstract 

We discuss the stringy properties of high-energy QCD using its hidden integrabihty in 
the Regge limit and on the hght-cone. It is shown that multi-colour QCD in the Regge limit 
belongs to the same universality class as superconformal M=2 SUSY YM with Nf = 2Nc 
at the strong coupling orbifold point. The analogy with integrable structure governing the 
low energy sector of M=2 SUSY gauge theories is used to develop the brane picture for 
the Regge limit. In this picture the scattering process is described by a single M2 brane 
wrapped around the spectral curve of the integrable spin chain and unifying hadrons and 
reggeized gluons involved in the process. New quasiclassical quantization conditions for the 
complex higher integrals of motion are suggested which are consistent with the 5— duality 
of the multi-reggeon spectrum. The derivation of the anomalous dimensions of the lowest 
twist operators is formulated in terms of the Riemann surfaces. 



1 Introduction 

During the recent years it was recognized that there is the hidden integrability behind the evo- 
lution equations in asymptotically free Yang-Mills theories |]1], 0. It was found that high-energy 
asymptotics of scattering amplitudes in QCD in the Regge limit are described by evolution equa- 
tions which are ultimately related to the SL{2, C) XXX Heisenberg integrable spin chains. More 
recently, a similar connection was observed for another physically interesting limit, namely the 
QCD dynamics on the light-cone. There, the logarithmic evolution of the composite light-cone 
operators is connected with the real SL{2,M.) XXX Heisenberg spin chain and the spectrum 
of anomalous dimensions of those operators was identified with the spectrum of this magnet 
[0, ^, ^ 1^. For the time being it seems that the phenomena is quite general and the effective 
QCD dynamics in different kinematical limits is integrable indeed. 
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On the other hand, there is another class of (supersymmetric) Yang-Mills theories in which 
integrability emerges. In the A/'=2 SUSY YM theory which can be solved in the low-energy 
limit exactly ^ the relevant integrable system was identified as the complexified periodic Toda 
chain. The low-energy effective action and the BPS spectrum in the theory can be described in 
terms of the classical Toda system [^]. When one adds the fundamental matter with Nf < 2Nc, 
which does not change the asymptotically free nature of the theory, the universality class of 
integrable model is changing. Instead of the classical Toda system one gets the classical XXX 
spin chain 0. In the case of the superconformal Nf = 2Nc theory, the solution leads to the 
periodic XXX chain in the magnetic field |]TU[. Because these magnets are classical, the spins are 
not quantized and their values are determined by the matter masses. In the special case, when 
the matter is massless, one arrives at the classical XXX magnet of spin zero. 

At the first glance there is nothing in common neither between high-energy (Regge) limit of 
QCD and low-energy limit of SUSY YM, nor between the corresponding integrable models. In 
this paper we shall demonstrate the opposite and show that there is a deep relation between the 
Regge limit of QCD and some particular limit of superconformal A/'=2 SUSY YM at Nf = 2Nc. 
This relation is based on the fact that in the both cases we are dealing with the same class of 
integrable models - the SL(2, C) homogenous spin zero Heisenberg magnets. However, there is 
the key difference between two theories. The integrable model in the SUSY YM is a classical 
Heisenberg magnet, whereas the Regge limit of QCD is described by a quantum Heisenberg 
magnet. To bring two pictures together, it is natural to develop quasiclassical description in 
QCD case and quantize the classical picture in SUSY YM. The first issue was considered in 
[PI , [T2| |13| while some arguments concerning the quantization of the integrable system related 
to SUSY YM theories can be found in 

The central role in our approach will be played by Riemann surfaces which appear naturally 
as spectral curves in a description of classical spin chains. To some extent the Riemann surface 
fixes the universality class of the classical integrable model and it defines general solution to 
the classical equations of motion. This solution is a finite-gap soliton expressed in terms of 
6'— functions on the Riemann surface. This surface has a natural complex structure whose moduli 
are parametrized by the integrals of motion. The classical dynamics is linearized on the Jacobian 
of this Riemann surface. In the case of the SUSY YM these Riemann surfaces encode the 
expression for the low energy effective actions |P, 1^. The moduli of the complex structure are 
fixed by the vacuum expectation values of the scalar in M=2 theory. One gets Riemann surfaces 
in the quasiclassical limit of the Schrodinger equation for the quantum spin chain related to 
QCD in the Regge limit |13[. The main goal of this paper is to use these Riemann surfaces to 
develop the stringy representation for the effective dynamics of QCD in the Regge limit exploring 
the known stringy realization of low-energy regime of SUSY YM. 

Our starting point is the observation that the Seiberg-Witten curve for the superconformal 
M=2 SUSY YM theory with gauge group SU{N^ and the number of flavours Nf = 2Nc coincides 
with the spectral curve for Heisenberg spin magnet describing multi-colour QCD in the Regge 
limit. On the SUSY YM side, we still have several unidentifled parameters. The flrst one is 
the rank of the gauge group SU{Nc). On the QCD side, in the multi-colour limit, we also have 
one additional parameter, namely the number of reggeized gluons N = 2, 3, ... exchanged in the 
channel [13, |TB|, 0. It will be shown later that the correspondence between these two theories 
implies that the number of reggeized gluons = A^^^. However the superconformal theory has 
additional freedom due to the arbitrary value of the bare coupling constant tq = + 2tt9 and 
the Seiberg-Witten curves are bundled over the complex tq plane. On the complex tq— plane 
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there are singular points at which the discriminant of the curve vanishes, i.e. the spectral curve 
becomes degenerate. It turns out that the Regge limit of the multi-colour QCD is related to one 
of these singular points, the so-called strong coupling orbifold point [|18|, |19| which is S"— dual to 
another singular point corresponding to the weak coupling e^'^*'^" —>■ regime. Let us emphasize 
that in the superconformal theory the effective coupling constant does not coincide with the bare 
one, because there are finite perturbative one-loop and nonperturbative instanton corrections to 
the low-energy effective action. At the strong coupling regime under consideration one has to 
take into account contributions from all numbers of instantons. 

After identification of the universality class one can investigate relevant S"— duality properties 
of both theories, which can be read off from the spectral curve. In the SUSY case, the S"— duality is 
the famous electric-magnetic duality. The Seiberg-Witten solution follows from the compatibility 
of the duality with the renormalization group (RG) flow of the effective action on the moduli 
space of the theory. In the QCD case the situation is more subtle since contrary to the SUSY 
case we have to deal with the quantum dynamics of the same integrable system and one has 
to formulate the precise meaning of the duality at the quantum level. It turns out that the 
duality manifests itself in the spectrum of the integrals of motion We formulate a new 

WKB quantization conditions obeying the duality property and demonstrate that their solutions 
are in a good agreement with the numerical calculations of the spectrum of the three reggeon 
system [^]. This strongly support the conjecture that the notion of duality survives at the 
quantum level. In result we predict that the quantum spectrum of the iV— reggeon system in the 
multi-colour QCD is parametrized by two {N — 2)— dimensional vectors n and fh, which have the 
meaning of "electric" and "magnetic" quantum numbers. 

We reformulate the above mentioned WKB quantization conditions in stringy terms using 
the brane realization of the BPS spectrum in superconformal M=2 SUSY YM theory pT| , 



Formally it is a simple task to identify the corresponding geometrical objects however the proper 
meaning of the "magnetic" quantum number in the Regge limit is a subtle point. The duality 
combined with RG behaviour give rise to brane realization of the ground state in the SUSY YM 
theories. There are two apparently different IIA/M and IIB/F pictures. In the IIA/M picture 
the ground state is described by the M5 brane with the world-volume i?^ x S where S is the 
spectral curve of the corresponding integrable system [|21], The stable BPS states correspond 
to the M2 branes ending on this M5 brane. In the IIB/F picture |22] one has to consider a D3 
brain in the orientifold background and BPS states are described by the dyonic strings stretched 
between the D3 brane and orientifold. 

In this paper we shall start to develop a similar brane picture for the Regge limit of the 
multi-colour QCD. We shall argue that the dynamics of the multi-reggeon compound states in 
multi-colour QCD is described in this picture by a single M2 brane wrapped around the spectral 
curve of the Heisenberg magnet. The fact that we are dealing with the quantum magnet implies 
that the moduli of the complex structure of the spectral curve can take only quantized values. 
We shall present a raw of arguments supporting this picture. Let us note that in the Regge limit 
there is a natural splitting of the four-momenta of reggeized gluons into two transverse and two 
longitudinal components, which is unambiguously determined by the kinematics of the scattering 
process. One can make a Fourier transform with respect to the transverse momenta and define 
a mixed representation: two-dimensional longitudinal momentum space and two-dimensional 
impact parameter space for which it is natural to use complex coordinates z and z = z* . The 
Riemann surface which the brane is wrapped around lives in this "mixed" coordinate-momenta 
space. Different projections of the wrapped brane represent hadrons and reggeons. The genus 
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of this surface is fixed by the number of reggeons. Let us note that muhi-reggeon states appear 
from the summation of only planar diagrams in QCD and from the point of view of topological 
expansion they all correspond to the cylinder-like diagram in a colour space. It contribution to 
the scattering amplitude is given by the sum of iV— reggeon exchanges in the t— channel with 
= 2, 3.... It is amusing that in result we get a picture where we sum over Riemann surfaces of 
arbitrary genus g = N — 2. 

Finally, we make some comments on a dual (super)gravity realization of the Regge limit of 
multi-colour QCD. Our approach is different from the recent studies p4| , |25| , [26[| of high-energy 
scattering in QCD based on the AdS^ geometry within the AdS/CFT correspondence E7^. We 



speculate that a reggeized gluon in QCD can be identified as a singleton in the AdS^ space. 

The paper is organized as follows. In Section 2 we remind how the integrable spin chains 
appear in the Regge limit of QCD. In Section 3 we review the main features of the finite-gap 
solutions to the classical equation of the motion in the periodic spin chains which are relevant 
for our consideration. We also discuss their quantization using the method of Separated Vari- 



ables [2S]. In Section 4 we identify the universality class of the Regge limit of multi-colour QCD 



by comparing the Riemann surfaces in QCD and A/'=2 SUSY YM. Based on this identification, 
we make a conjecture about brane realization of QCD in the Regge limit. In Section 5 we investi- 
gate the duality properties of the multi-reggeon states and provide WKB like description of their 
spectrum in terms of Riemann surfaces. In Section 6 we propose the stringy/brane picture for 
the description of the anomalous dimensions of the lowest twist operators in QCD. In Section 7 
several comments concerning the supergravity dual description of the Regge limit are presented. 
The obtained results are discussed in the Conclusion. 



2 Regge asymptotics in QCD and Heisenberg magnets 

It is well known from the old days of the Regge theory that the scattering amplitudes of hadrons 
grow at high-energy. Useful framework for understanding this phenomenon in QCD comes from 
the studies of the asymptotic behaviour of the scattering amplitude of two hadronic states con- 
sisting of a pair of heavy quarks |2^ , the so-called onia system. The advantage of such system is 
that it captures all essential features of the Regge phenomenon in QCD and, at the same time, 
the scattering amplitudes can be calculated perturbatively order-by-order in the QCD coupling 
expansion. Denoting the invariant mass of the scattering onia by Q"^ and their center-of-mass 
energy by s, we introduce the dimensionless Bjorken scaling variable x = Q'^/s. The Regge limit 
corresponds to the asymptotics at large s = Q'^/x with fixed Q^, or equivalently to small x. 

One expects that the calculation of the high-energy asymptotics of the onium-onium cross- 
section should eventually match the Regge model prediction. The Regge model interprets the 
growth of the scattering amplitudes at high energy s, or equivalently at small x, by introducing 
the notion of the Pomeron as the Regge pole with the quantum numbers of the QCD vacuum. 
Its contribution of the onium-onium cross-section at small x takes the form 

aUx,Q')=x'^''P-'^(3f{Q')(3f{Q'), (1) 

with a-p the Pomeron intercept and P^siQ'^) residue factors corresponding to the onia A and 
B, respectively. Despite the fact that the Regge model provides a successful phenomenological 
description of the experimental data in terms of "hard" (perturbative) and "soft" (nonperturba- 
tive) Pomerons its status within QCD remains unclear. 
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The first attempts to understand the "hard" Pomeron within perturbative QCD have been 
undertaken more than 20 years ago and they have led to the discovery of the BFKL Pomeron ||30|1 . 
In the BFKL approach, two onia scatter each other by exchanging soft gluons in the t— channel. 
Calculating the corresponding Feynman diagrams in powers of as one obtains the following 



general form of the perturbative expansion of the cross-section at small x and fixed Q , |15, 16, 1 / 



oo 

ctUx,Q'') = aBorn J](asiVc)'"-"(asiVc In x)"/„,„(Q^) , (2) 



m=0 n=0 



with o"Born being the cross-section at the Born level. Here, the double sum involves two parameters 
of the expansion - the QCD coupling constant, ctg <^ 1, and the energy dependent parameter 
as In 1/x that becomes large at small x. The coefficient functions fm,n{Q'^) depend on the invariant 
mass of the onia, Q^, as well as on the number of colours, l/N^- They can be calculated 
perturbatively for large invariant mass ^ ^qcd making use of the wave functions of the 
onia states. The number of different terms in rapidly increases to higher orders in and in 
order to find atot{x, Q"^) at small x one has to develop a meaningful approximation which, firstly, 
correctly describes the small— x asymptotics of the infinite series @ and, secondly, preserves the 
unitarity constraint 

crtot(2;,<5^) < const. X ln^(l/a;) (3) 

as X — > 0. To satisfy the first condition, one can neglect in (Q) the terms containing fm,n with 
n < m as suppressed by powers of as with respect to the leading term fm,m- The resulting series 
defines the onium-onium cross-section in the leading logarithmic approximation (LLA). It can 
be resummed to all orders in as leading to the BFKL pomeron |3( 



°° ™-asAfc41n2/7r 

m=o ^/asNc\nl/x 

The BFKL Pomeron leads to unrestricted rise of the cross-section and violates the unitary bound 
(^. In order to preserve the unitarity of the S-matrix of QCD and fulfill (^), one has to take 
into account an infinite number of nonleading terms in (0). This means that with the unitarity 
condition taken into account the series does not have small parameter of expansion like as- 
Instead of searching for this parameter one may start with the LLA result and try to identify 
the nonleading terms in (H), which should be added to (§) in order to restore the unitarity. In this 



way, one arrives at the generalized leading logarithmic approximation [|T^, |T5|, |T6[. One should 
keep in mind, however, that in this approximation the unitarity is preserved only in the direct 
channels of the process but not in subchannels corresponding to the different groups of particles 
in the final state. 



2.1 Generalized LLA 

In the generalized LLA, one uses the remarkable property of gluon reggeization in order to 
formulate a new diagram technique for calculation of the onium-onium scattering amplitudes 
[0, [1^ . Namely, an infinite set of standard Feynman diagrams involving "bare" gluons can 
be replaced by a few reggeon diagrams describing propagation of reggeized gluons, or reggeons, 
and their interaction with each other. Each reggeon diagram appears as a result of resummation 
of an infinite number of Feynman diagrams with "bare" gluons. In this way, one defines an 
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effective tlieory in wfiicfi tfie reggeized gluons start play a role of a new elementary fields. In this 
theory, the scattering amplitudes describe the effects of the propagation of the reggeized gluons 
in the t— channel and their interaction with each other. As usual, each diagram in the effective 
theory is equivalent to an infinite sum of the original QCD Feynman diagrams. 

The resulting set of the reggeon diagrams defining the onium-onium cross-section in the 
generalized LLA is shown in Figure 1. According to the optical theorem, the cross-section is 
given by an imaginary part of their contribution. The vertical curly lines represent the reggeons 
propagating in the t— channel and the triple gluon vertices describe the effective BFKL interaction 
between the reggeons and the gluons propagating in the s— channel. Upper and lower blobs denote 
the coupling of the reggeons to two onia states A and B. 




Figure 1: The Feynman diagrams contributing to the scattering amplitude of the onium-onium 
scattering in the generalized LLA and describing the reggeon exchange in the t— channel. 



One of the peculiar features of the reggeon scattering in the generalized LLA is that it is 
elastic and pair-wise. The number of reggeons exchanged in the t— channel is not changed. 
This allows us to interpret the diagrams shown in Figure 1 as describing the propagation of the 
conserved number N = 2, 3,... of pair-wisely interacting reggeons. At = 2 the diagram in 
Figure 1 determines the leading logarithmic result for the scattering amplitude, Eq. (^, and it 
contributes to the coefficient functions /m,m(Q^) in (H) to all orders of perturbation theory. The 
contribution of the diagram with A^ = 3 is suppressed by a power of as with respect to that for 
N = 2 and it determines the first nonleading coefficient fm,m-i in (0). 

The calculation of the A^— reggeon diagrams can be performed using the Bartels-Kwiecinski- 
Praszalowicz approach [^, |T^ . The diagrams shown in Figure 1 have a generalized ladder form 



and for fixed number of reggeons, A^, their contribution satisfies Bethe-Salpeter like equations 
for the scattering of A^ particles. The solutions to these equations define the colour-singlet 
compound states {"^n) built from A^ reggeons. These states propagate in the t— channel between 
two scattered onia states and give rise to the following (Regge like) high energy asymptotics of 
the the onium-onium cross-section 

. „—asNcei^ /-IT 

aUs)- E (a^Nr ^—^-^ Pf\Q')P'B\Q') (5) 
Af=2^,... y/asNcinl/x 

with the exponents ejy defined below. Here, each term in the sum is associated with the con- 
tribution of the A^— reggeon compound states. At A^ = 2, the corresponding state defines the 
BFKL pomeron M) with 62 = 41n2. 
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The TV— reggeon states are defined in QCD as solutions to tlie Schrodinger equation 

(y. 

7iiV^(£l, Z2, Zn) = —NcSN'^izi, Z2, Zn) (6) 

TT 

with the effective QCD Hamiltonian Ti^ acting on two-dimensional transverse coordinates of 
reggeons, {k = 1, ...,N) and their colour SU{Nc) charges 

^jv = -^ J2 Hijt^t"-. (7) 

l<i<j<Af 

Here, the sum goes over all pairs of reggeons. Each term in this sum is factorized into the product 
of two operators acting on the color indices and the two-dimensional coordinates. The former is 
given by the direct product of the gauge group generators in the adjoint representation of the 
SU{Nc) group, acting in the color space of z— th and j— th reggeons, 

= I®...(g)e ^...(^I, {t'')bc = -ifabc, 

i 

with fabc the structure constants of the SU{Nc). The operator Hij acts only on the transverse 
coordinates of reggeons, Zi and zj. The Hamiltonian commutes with the total color charge 
of the system of iV— reggeons. Solving (H), we are looking for the colour- singlet states, that is 
the eigenstates with the total color charge equal to zero, X^ili ''-i'l^Af) — 0- The residue factors 
in d^) measure the overlap of Z2, zn) with the wave functions of two scattered particles, 
P^^\Q^) = {A\^n) and similar for 

To get some insight into the properties of the A^— reggeon states, it proves convenient to 
interpret the Feynman diagrams shown in Figure 1 as describing a quantum-mechanical evolution 
of the system of particles in the t— channel between two onia states \A) and \B) 

atot(s) = 5^(asiVe)^(A| ei°(V-)-^- \B) , (8) 

N>2 

with the rapidity Inx = In Q^/s serving as Euclidean evolution time. To find the high energy 
asymptotics of , one has to solve the Schrodinger equation (^) for arbitrary number of reggeized 
gluons N, expand the operator over the complete set of the eigenstates of Hn and, then, 

resum their contribution for arbitrary A^. In this way, one arrives at (^. 

It follows from that for fixed number of the reggeons, N, and large evolution time, 
Inl/x ^ 1, the dominant contribution to atotls) comes from the eigenstates of Ti^ with the 
maximal energy. This allows us to identify the exponent (— ^tv) in (^) as the energy of the 
ground state of the Hamiltonian {—Hn)- It turns out that the energy spectrum of the reggeon 
Hamiltonian is gapless and the contribution of the eigenstates next to the ground state amounts 
to the appearance of the additional factor ~ (og In l/x)"^/^ in the r.h.s. of (|^). The contribution 
of the A^— reggeon state to crtot(s), Eq. (^, depends crucially on the sign of the energy - for 
e^v > (or e^r < 0) it increases (or decreases) as s — oo. 

Finding the spectrum of the compound states (^ for arbitrary number of reggeized gluons A^ 
and eventual resummation of their contribution to the scattering amplitude (^ is the outstanding 
theoretical problem in high-energy QCD. At A^ = 2, the solution to (P) has been found a long 
time ago - the well-known BFKL Pomeron |Q, £2 = 4 In 2. At A^ = 3 the solution to (|^) - 



the Odderon state in QCD |31||, was formulated only a few years ago £3 = —0.2472.., by 
making use of the remarkable properties of integrability of the effective QCD Hamiltonian [||, 0. 
Recently, a significant progress has been achieved in solving (0) for arbitrary number of reggeons 



in the milti-colour limit |3^, |3^, ^ . 
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2.2 Large A^^ limit 

At large Nc one can expand the matrix elements in the r.h.s. of @ in powers of l/N"^ and obtain 
a topological expansion of crtot(s) "^^^ leading term of this expansion corresponds to the 

planar cylinder diagram in which the top and the bottom disks correspond to two onia states 
and the walls are formed by N reggeons. In this diagram, in order to preserve the planarity, each 
reggeon is allowed to interact only with two neighbouring reggeons. Going back from the planar 
diagrams to the Hamiltonian TYat, this leads to simplification of the reggeon interaction in the 
multi-colour limit. Namely, for a^N^ = fixed and ^ oo, the Hamiltonian (|^) can be written 
as 

= ^ E ^^.^+1 + ^(^c') (9) 

l<i<N 

with H^^j^+i = Hjy i. Thus, the pair- wise interaction between reggeons in is reduced in the 
multi-colour limit to a nearest-neighbour interaction and, in addition, the colour operator is 
replaced by a c— valued factor —Nc/2. As was already mentioned, the two-particle reggeon 

Hamiltonian -ffj.i+i acts on the two-dimensional (transverse) reggeon coordinates z = {x, y) and 



it coincides with the BFKL kernel ||3^. It becomes convenient to introduce the complex valued 
(anti)holomorphic coordinates, z = x + iy and z = x — iy, and parameterize the position of the 
kth reggeon as Zk = {zk,Zk). One of the remarkable features of the BFKL kernel is that it is given 
by the sum of two mutually commuting operators acting separately on z— and coordinates 
of the reggeons. Making use of this property, one can rewrite the iV— reggeon Hamiltonian (j^) 
as 
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= ^^Hr, + Hn) + 0{N~') . (10) 

Here the hamiltonians and Hn act on the (anti)holomorphic coordinates and describe the 
nearest-neighbour interaction between reggeons 

N N 
m=l m=l 

with the periodic boundary conditions z^+i = Z\ and = Z\. The interaction Hamilonian 
between two reggeons with the coordinates [z^^ ^i) and (^2, ^2) in the impact parameter space, is 
given by the BFKL kernel 

i7(zi, Z2) = -^( J12) - ^(1 - J12) + 2V^(1) , (11) 
where ^/'(x) = d\nr{x)/dx and the operator J12 is defined as a solution to the equation 

Jl2{Jl2-l) = -{^l-Z2yd,d2 (12) 

with dm = d/dzm- The expression for H{zi,Z2) is obtained from (|TTp by substituting Zk Zk. 

Since the Hamiltonian and Hn act along two different directions on the two-dimensional 
plane, they commute with each other, [Hn, Hn] = 0. This allows us to reduce the original 
Schrodinger equation (^ to the system of two one- dimensional Schrodinger equations for the 
Hamiltonians Hn and Hn- For fixed number of reggeons, A^, the operator Hn describes the 
system of A^ particles on a complex 2— line interacting with their neighbours through the two- 
particle hamiltonian (|lT]). This quantum- mechanical system can be solved exactly for the system 
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with N = 2 particles leading to the BFKL Pomeron. It is not obvious, however, whether the 
exact solution can be found for an arbitrary number of reggeons. For this to occur, the system 
of N reggeons should possess an additional symmetry. It turns out that such a hidden symmetry 
indeed exists [|l], @] . Namely, the Schrodinger equation for the A^— reggeon system in multi-colour 
limit contains the set of the mutually commuting integrals of motion and {k = 2,...,N) 
defined as 

Qk = i'^ ^ Zjij2---^jk-i,jk^jk,ji^Zj-^---dzj^_^dzj^ (13) 

^<jl<j2<-<jk<N 

with Zjk = Zj — Zk- The charges qt are given by similar expressions in the sector. The 
eigenstates ^{zi, Z2, zn) have to diagonalize these operators and the corresponding eigenvalues 
q = {q^, qi- = with k = 2, N, form the complete set of quantum numbers parameterizing 
the spectrum of the Schrodinger equation (H), ejy = e]\f{q). This implies that the Schrodinger 
equation (P) for the A^— reggeon states is completely integrable in the multi-colour limit and it 
can be solved by applying the powerful Quantum Inverse Scattering Method well-known in the 
theory of integrable models. 

The famous example of the exactly solvable many-body quantum mechanical system is the 
one-dimensional XXX Heisenberg chain of N interacting spins described by the hamiltonian 



N 

rrS=l/2 _ ST- q q 



m=l 

where Sm are the SU{2) generators of the spin s = 1/2 acting in the m— th site and Stv+i = 5*1. It 
turns out p|, ^ that this simple model admits nontrivial generalizations for an arbitrary complex 
value of the spin s. In the latter case, the spin operators 5"^ are the generators of an infinite- 
dimensional representation of the SL{2, C) group of the principal series. Remarkably enough, 
in the special case of the SL{2, C) spin s = 0, the Hamiltonian of such completely integrable 
noncompact XXX Heisenberg magnet becomes identical to the QCD Hamiltonian describing 
the interaction between A^— reggeons in the multi-color limit. 

To explain this correspondence let us consider the one-dimensional lattice with periodic 
boundary conditions and with the number of sites, N, equal to the number of Reggeons. Each 
site is parameterized by the two-dimensional vector Zm = {zm, Zm) with Zm and z^ being holo- 
morphic and antiholomorphic coordinates, respectively, and m = 1, .., A^. In addition, one assign 
to each site the spins operators Sp^ and Sp^, which are the six generators of the principal series 
of the SL{2, C) group. They are realized on the quantum space of the model as the differential 
operators 

S', = zA^ + s, S,= -d., , = zld.^ + 2sz, . (14) 

The operators Sp^ are given by similar expressions with zt, s replaced by Zk, s = 1 — s*, 
respectively. The pair of complex parameters (s, s) specifies the SL{2, C) representation. For 
the principal series of the SL{2, C) they are parameterized by integer rig and real Ug as s = 
(1 -|- ns)/2 + iug and s = 1 — s*. To match the reggeon Hamiltonian, Eq. (P), the spins have to 
be equal to s = and s = 1. In this representation, the BFKL kernel defining the interaction 
between two reggeons describes the interaction between the nearest spins. The Hamiltonian of 
exactly solvable noncompact SL{2, C) Heisenberg magnet of spin (s, s) is defined as p|, ^ 



d 



Hjq — \ ^ Hm,m+1 ; -f^m,m+l — ~^ l^- -Rm,m+1 (^) ^) (15) 

m=l 
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where the operator Rm,m+i, the so-called fundamental i?— matrix for the SL{2, C) group, satisfies 
the Yang-Baxter equation and it is given by 



_ r(m)r(i + m) r(i - J12 -m)r(Ji2 - w) 

i?12(M,M) = . — X — • x-p/ T (16) 

r(— m)r(i — lu) r(i — Ju + lu) r(Ji2 + m) 

with J12 defined as the sum of two spins 

Jk,k+i{JKk+i - 1) = (5(^) + S^'^^'^f (17) 

with S'i^'^^'' = Sa \ and JA;,fe+i is defined similarly. The operator Ri2{u,u) acts on the tensor 
product V 1^ V with V = being the representation space of the principal series of the 

SL{2X)- 

Substituting (|1^) into ([T5|), we verify that the Hamiltonian of noncompact Heisenberg magnet 



of the SL{2, C) spin s = and s = 1 coincides with the Hamiltonian of Reggeon interaction ([Tl|) . 
This leads to the following identity between two Hamiltonians 

47r 

Thus, the problem of finding the spectum of the A^— reggeon states in multi-colour QCD is 
reduced to solving the Schrodinger equation for a peculiar completely integrable model - the 
noncompact Heisenberg magnet with spin operators being the generators of the principal series 
of the SL{2, C) group of the spin s = 0. This model has a number of remarkable properties 
[@, which make it completely different from conventional Heisenberg spin magnets studied 
so far. Firstly, the SL{2, C) representation of the principal series does not have the highest 
weight and, as a consequence, the Algebraic Bethe Ansatz is not applicable to construction the 



eigenstates of the model. Secondly, despite the fact that the Hamiltonian (|T0|) can be split 
into the sum of holomorphic and antiholomorphic mutually Hamiltonians acting on the z— and 
coordinates of the particles, the two sectors are not independent. The interaction between 
two sectors occurs through the condition for the wave function to be a single valued function 
on the z= (2;, ^)— plane. This requirement plays a crucial role in establishing the quantization 
conditions on the integrals of motion and in finding the energy spectrum of the model p 



3 Finite gap solutions and their quantization 

To get some insight into the properties of the Schrodinger equation (|^), it becomes convenient to 
consider a classical analog of the noncompact Heisenberg spin magnet. As was shown in |]rT|, |T^, 
the corresponding classical Hamiltonian model naturally appears in the WKB solutions of the 
Schrodinger equation (|^). From point of view of classical dynamics, the system describes two 
copies of one-dimensional spin chains "living" on the complex z— and z— lines. Each of these 
spin chains can be analyzed separately and, for simplicity, we shall concentrate in this Section 
on the holomorphic system. 
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3.1 Classical dynamics 

Following the Quantum Inverse Scattering Method, the classical homogenous spin s = chain of 
the length N is defined in the holomorphic sector by the 2x2 Lax matrices 

Lkix)=x-l + ySt-a''=(''~. '2'^' 1 (19) 

where a" are the Pauli matrices and x is the spectral parameter. Here, S'l stand for classical spins. 
They are obtained from the spin operators in the holomorphic sector, Eq. (p!^), by substituting 
the derivatives with respect to 2;— coordinates by the corresponding classical momenta idz^ Pk 



Qk 



classical 

l<jl<j2<---<Jfe<Af 



In this way, one can define on the phase space of the model the Poisson brackets 

{zk,Pj} = 6jk (21) 

and verify that the Poisson brackets of the dynamical variables Sa, a = 1,2,3 (taking values in 
the algebra of functions) look like 

{S^,,S'^} = -ze,,,S^,6,k, (22) 

so that {Sa} plays the role of angular momentum ("classical spin") giving the name "spin-chains" 
to the whole class of systems. Algebra ( p^ ) has an obvious Casimir function = X]a=i ^j^j ' 
which does not depend on the number of the site, j = 1, ...,N for the homogenous spin chain. 
In similar manner, one can define the classical analogs of the operators qk and the Hamiltonian 
Hn- 

The equations of motion of the classical spin magnet defined in this way look like 
They possess the set of — 1 conserved charges g„ 

dtQn = {qn,nN} = . (24) 

with n = 2,...,N and their solutions define the reggeon trajectories Zk = Zk(t) subject to the 
periodicity condition Zk+N(t) = Zkif) and pk+Nif) = Pk(t). 

To demonstrate a complete integrability of the evolution equations, one observes that the 
system (^) is equivalent to the matrix Lax pair relation 

dtLk{x) = {Lk{x),nN} = Ak+i{x)Lk{x) - Lkix)Akix) , (25) 

where Ak{x) is a 2 x 2 matrix depending on the momenta and the coordinates of the particles. 
Let us introduce the Baker- Akhiezer function \l//c(x;t) as a solution of the following system of 
matrix relations 

Lkix)^kix,t) = ^k+iix,t), dt'^kix^t) = Akix)'^kix,t) . (26) 
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The two-component Baker-Akhiezer function defined in this way depends on the site number k. 
The Lax operator Lk{x) shifts it into the neighbouring site. One can introduce the monodromy 
operator producing the shift of along the whole chain as 

^k+N{x;t)=TN{x)^k{x;t), Tn{x) = Ln{x) . . . Li{x) (27) 

The periodic boundary conditions on the solutions to (|23|) can be easily formulated in terms of 
the Baker-Akhiezer function as 

^k+Nix;t) =w^kix;t) (28) 

where w is the Bloch-Floquet factor. According to (^) and (^), w is the eigenvalue of the 
monodromy operator Tjv(x), and therefore it satisfies the characteristic equation 

det(T7v(x) -wl)=0 (29) 

From Eqs. (|25|) and (^), the time dependence of the monodromy operator is given by dtTj\i{x) = 
y4i(x)T/v(x) — T/v(x)y4i(x), and, as a consequence, its eigenvalues w are conserved quantities 
generating a complete set of integrals of motion. Evaluating the l.h.s. of (^), one obtains the 
spectral curve equation 

-w -TrTNix) + detT{x) = 0. (30) 
Using the definitions (^) and (|19[) one gets 

tN{u) = TiTn{x) = 2x^ + q2X^-^ + ... + q^-iX + (31) 

with qk being the integrals of motion, Eq. (|2Q|) , and 

N 

detT(x) = JJdetLfc(x) = x^^. (32) 

k=l 

Introducing the complex function y{x) = w — x"^^ /w one obtains the equation of the spectral 
curve in the form 

: y'' = tUx) - 4x2^ = (4x^ + gs^^"' + ... + g^)(g2X^-' + ... + Qn) ■ (33) 

For any complex x in general position the equation ( |30D has two solutions for w, or equivalently 
yix) in dH). 

Under appropriate boundary conditions on \l/fe(a;) these solutions define two branches of the 
Baker-Akhiezer function, \l/^(x). Then, being a double- valued function on the complex x— plane, 
\l/fc(a;) becomes a single- valued function on the Riemann surface corresponding to the complex 
curve Fat. This surface is constructed by gluing together two copies of the complex x— plane 
along the cuts [62,63], [e2N-2,ei] running between the branching points ej of the curve (|33|) , 
defined as simple roots of the equation 

4(6,)=46f, j = l,...,2Ar-2. (34) 

Here positions on the complex plane depend on the values of the integrals of motion q2, qs, 
q^. In general, the Riemann surface defined in this way has a genus g = N — 2 which depends 
on the number of reggeons, N. For instance, it is a sphere at = 2 and a torus at = 3. 
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Applying the standard methods of the finite-gap theory, it becomes possible to construct the 
explicit expression for the Baker- Akhiezer function in terms of the theta-functions defined on 
the Riemann surface (|33| ) and, finally, obtain the solutions to the classical equations of motion 



3[). The explicit expressions can be found in [12|. So far we restricted our consideration to the 



classical dynamics in the z— sector. Obviously, similar relations hold in the 2— sector. Since the 
integrals of motion, and qk, and the Hamiltonians in two sectors, and Hn, are conjugated 
to each other, the solutions of the classical equations of motion in two sectors can be obtained 
one from another. 

Combining together the classical motion along the z— and ^—directions, we find that the 
resulting expressions describe the finite-gap soliton wave propagating in the system of particles 
located on the two-dimensional plane. In terms of the spectral curve the same classical 

djTiamics corresponds to the motion of particles on the Riemann surface (|33D. The motion on 
the Riemann surface is not confined to the kinematically allowed bands corresponding to the 
a— cycles and the classical trajectories wrap along both the a— and /9— cycles. This is one of the 
manifestation of the fact that we are dealing with two-dimensional system. As we see below, it 
will play an important role in our analysis. 

3.2 Quantum case 

Let us now consider the quantum evolution of the system with the Hamiltonian defined in ([T5|). 
As was explained in the Section 2.1, the eigenstates of the noncompact SL{2,C) Heisenberg 
chain of length N and spin (s = 0, s = 1) play the role of the wave functions of the colourless 
compound states of A^ reggeized gluons. In addition, the ground state of the system controls the 
high-energy asymptotics of the scattering amplitudes, Eq. (H). 

Due to complete integrability of the system, the A^— Reggeon states have to diagonalize si- 
multaneously the integrals of motion and q^ with k = 2, ...,N Their corresponding eigenvalues 
provide the total set of the quantum numbers of the A^— Reggeon states, (g, g). In particular, 
the "lowest" integrals of motion, q2 and q2, coincide with the quadratic Casimir operators for 
the total SL{2, C) spin of the system and their eigenvalues are related to the conformal SL{2, C) 
spin of the A^— Reggeon state as 

q2 = -h{h - 1) , q2 = ~h{h - 1) , q2 = q*2, (35) 

where all possible values of h and h = 1 — h* can be parameterized by integer n and real u 

h = —^ n = Z, z/ = lR. (36) 

As to remaining charges, gs, qN and % = <iki their possible values are also quantized. The ex- 
plicit form of the corresponding quantization conditions is more complicated and it was obtained 



m 



According to (13), the integrals of motion are given by the A;— order differential operators (with 
k = 2,...,N) acting on the holomorphic and antiholomorphic reggeon coordinates. Therefore, 
their spectral problems leads to a complicated system of coupled differential equations on the 
wave functions of the A^— reggeon states \E'g^g(zi, z^)- The solution to this system can be found 
by applying the method of separated variables developed by Sklyanin [^. Namely, there exists 



the unitary transformation Ujs,x{z) 

t^pA,...,xjv_i(^i, •••,^7v) = (il, •••,^7V-l) , (37) 
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that allows us to transform the wave function from the original representation defined by the 
N coordinates of Reggeons on the two-dimensional plane to the representation of the separated 
coordinates x = (xi, ...,XAr_i) and p being the total momentum of the system. The eigenstate of 
the Hamiltonian in this representation is given by 

d'^Z f/p',xi,...,a;iv-i(^) , (38) 



where d'^z = (Pzi...(Pzn- The remarkable feature of the unitary transformation U^^g{z) is that 
it transforms the wave function q}(z) satisfying the original Schrodinger equation (P) into 
the wave function in the separated coordinates, ${g g}(x), which is factorized into the product of 
Q— functions depending on a single separated coordinate Xk and satisfying the Baxter equation 
(see Eqs. (|1]) and (g) below) 

^{q,q}{Xi, ...,XN-i) = Q(xi, Xi) . . . Q{xn-1,Xn^i) ■ (39) 

Here, the notation was introduced for the holomorphic and antiholomorphic components of the 
separated coordinates Xk = {xk,Xk)- The explicit construction of the transformation to the 
separated coordinates, Eq. (p8|), can be found in [Q. Requiring U^',cci,...,xn_i{z) to be a single- 
valued function on the i*— plane, one finds that the possible values of the separated coordinates 
Xk have the following form 

ink _ , ink 

Xk = ^k Xk = Uk + — (40) 

with z/fe real and nk integer. 

The functions Q{x, x) entering (|39|) have to satisfy the functional Baxter equations in the 
x— sector 

x^Q{x + i,x) + x^Q{x — i,x) = tN{x)Q{x, x) (41) 
and similar relation holds in the x— sector 

{x + i)^Q{x + i,x) + {x- i)^Q{x -i,x) = tN{x)Q{x, x) . (42) 

Here, t]sf{x) is a polynomial of degree iV in x with the coefficients given by the eigenvalues of the 
integrals of motion, Eq. (|3TD. tiv(x) is given by similar expression with qk replaced by qk = oX- 
To make the correspondence with the classical mechanics, one introduces the shift operator on 
the space of functions Q(x, x) 

e^PQ{x,x) = Q{x±i,x), e^^ Q{x,x) = Q{x,x ±i) . (43) 

Then, the Baxter equations can be rewritten as 

[x^(eP + e-P) -t^,(x)] Q{x,x) = [x^(eP" + e-P") - t^(x)] x^Q(x,x) = 0. (44) 

Based on the Baxter equations, it is now possible to establish the correspondence between the 
quantum and the classical systems. Let us forget for a moment that p and x are operators 
in (^) and treat them as c— numbers. If one identifies the Bloch-Floquet factor in (|30| ) as 



w{x) = X exp(p), then it is easy to see that the l.h.s. of (^if) coincides with the equation of 
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spectral curve (|30|) . Going back from the classical to quantum system, one concludes that the 
Baxter equation (^4|) arises in result of the "quantization" of the spectral curve (|30|) on the space 
of the function Q{x,x). In other words, the constraint imposed by the spectral curve on the 
phase space of the classical system becomes an operator equation on the wave function of the 
same system after quantization. Thus the Riemann surface naturally enters into the quantum 
spectral problem. The important difference with respect to the classical case is that now the 
integrals of motion defining the moduli of the complex structure on (^) must be quantized. 

The Baxter equations do not fix the function Q{x,x) uniquely as their solutions are defined 
up to multiphcation by an arbitrary periodic function f{x + i,x) = f{x,x + i) = f{x,x). To 
avoid this ambiguity one has to impose additional conditions on the solutions to (|41|) and (|42|). 
As was demonstrated in |2D|, these conditions follow from the explicit expression for the 
kernel of the unitary transformation (^) and they can be formulated as requirement for Q{x,x) 
to have correct analytical properties and asymptotic behaviour at infinity. Namely, choosing 
X = X — in/2 and x = X + in/2 in accordance with (POj), and extending the values of A from real 
axis to the complex plane, one requires that Qq^g{x, x) should be a meromorphic function of A for 
fixed integer n with an infinite set of poles of the order not higher than N situated at the points 

{x^ = =Fim , = =F2(m - 1)} , m,m = 1,2,... (45) 

with m and rh being positive integer. In addition, the asymptotic behaviour of the solutions to 
the Baxter equation at infinity is given by 

g(A - m/2, A + m/2) ^-"^ e*® A^+^-^ + e"*® A^-^+i"^"^ (46) 
with the phase 6 = 0/i('?,^) depending on the quantum numbers of the state and the total 



SL{2, C) spins h and h = 1 — h* defined in (|35|). Up to an overall normalization, the above two 



conditions fix uniquely the solutions to the Baxter equation, and, in addition, they allow us to 
determine the spectrum of quantized integrals of motion and qk- 

Using the solution to the Baxter equations, we construct the wave function of the A^— reggeon 
states in the separated coordinates (|5^). To find the corresponding value of the energy one has to 



apply the resulting expression for the wave function to the Hamiltonian Hn transformed to the 
separated coordinates. This leads to the following remarkably simple expression for the energy 
of the iV— reggeon state 

Hn = i^\r^{x'^^ [Qix-i,x)]*Q{x + i,x)) . (47) 
dx x=o 

The relations (^) , (|4^) , ( ^5|) , (^) and (|47|) provide the basis for calculating the spectrum of the 
A^— reggeon states. 

The Baxter equations can be solved exactly at A^ = 2 and their solution is expressed in terms 
of the 3F2— hypergeometric series. Going though the calculation of the energy we arrive at 



E2[h) = 62 = -2 Re 

TT TT 



(48) 



This relation coincides with the well-known expression |^6[ for the energy of the N = 2 Reggeon 
compound state. The maximum value of the energy 

Er" = ^^41n2 (49) 

TT 
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defines the intercept of tlie BFKL Pomeron. For = 3 tlie solution to tlie Scfirodinger equation 
lias been first found in by direct diagonalization of the integral of motion q-^ in the original 
2— coordinates. Recently, the same result was obtained within the Baxter equation approach in 
[0, The ground state of the = 3 system defines the Odderon state in QCD |Q and its 
energy is given by 

"^^<=.e3 = _^. 0.247170... (50) 



771 max 
^3 



TT 



TT 



Contrary to (^91), the energy of the = 3 state is negative. This implies that its contribution 
to the scattering amplitude @ decreases at high-energy as s — cx), or x — 0. Therefore, it 
becomes of most interest to find the solution of the Baxter equations, Eqs. (^Tj) and (^), for 
higher > 4 Reggeon states and calculate the spectrum of the energies governing the high- 
energy asymptotics in (^. This problem has been recently solved in |l33| , pO| and the results for 



£7V are shown in the Figure 2. The detailed description of the spectrum can be found in 
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Figure 2: The dependence of the energy of the A^— reggeon states, E'^^^ = asNcSN/iT, on the 
number of particles A^. The exact values of the energy are denoted by crosses. The upper and 
the lower dashed curves stand for the functions 1.8402/(A^ - 1.3143) and -2.0594/(A^ - 1.0877), 
respectively. 



It was found in ||2^ that the A^— reggeon states have different properties for even and odd A^. 
For even A^, the ground state of the system is unique, whereas for odd A^ it is double degenerate. 
In the latter case, the color-singlet compound states built odd number of reggeized gluons have 
the same energy but different quantum numbers with respect to the charge conjugation. For odd 
A^ the A^— reggeon states with the charge parity C = 1 and C = —1 belong to the Pomeron and 
the Odderon sector, respectively. For even A^, the reggeon states have the charge parity C = 1. 

For odd N the ground state energy, e^, is negative and it increases with A^ approaching the 
value £2oo+i = from below. For even N, the ground state energy is positive and it decreases 
with A^ approaching the same value £200 = from above. We recall that for > (^at < 0) the 
contribution of the A^— reggeon state to the cross-section (^) increases (decreases) at high energy 
s. Thus, in the Pomeron and the Odderon sectors, the contribution of the A^— reggeon states to 
the cross-section ceases to depend on the energy s as A^ 00. It is interesting to notice that 
this result has been anticipated a long time ago within the bootstrap approach p5 |. 

As we have seen in the previous Section, the classical analogs of the A^— reggeon states are 
described by the finite-gap soliton waves propagating in the system of A^ particles on the plane. 
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The solitons are uniquely specified by the Riemann surface (^) whose moduh are defined by 

the conserved charges qk- The solution to the Schrodinger equation for the A^— reggeon states 
leads to the quantization of these charges. This allows us to interpret the whole quantization 
procedure described in this Section as quantization of the moduli space of the Riemann surface 



4 Universality class of the Regge limit 

In the previous Section the Riemann surfaces appeared within the context of high-energy QCD as 
auxiliary objects which were introduced to solve the equations of motion of the classical system 
of reggeized gluons and the Baxter equations for the corresponding quantum system. In this 
Section we shall argue that these surfaces have a definite physical meaning. Actually the situation 
is parallel to SUSY YM case where the auxiliary Riemann surface later on was identified as a part 
of the six-dimensional worldvolume of the M5 brane. The remaining four dimensions provide the 
worldvolume of the theory under consideration. 

Here we shall develop a similar brane picture for the Regge limit of multi-colour QCD. To do 
this we explore the analogy with the J\f=2 SUSY YM where the stringy/brane picture naturally 
emerges from the hidden integrability governing the low-energy effective actions. We shall argue 
that in the Regge limit the Riemann surface which has been considered earlier plays a similar 
role providing the worldvolume to the M2 brane describing the scattering process. In this way 
this Riemann surface fixes the universality class of the multi-colour QCD in the Regge limit. 

4.1 Riemann surfaces and QCD versus SUSY YM 

Let us recall the main features of the low-energy effective action in the M=2 SUSY YM theories 
relevant for our purposes [0, |3^. The key point is that the theory has the nontrivial vacuum 
manifold since the potential involves the term 



Here is the complex scalar field which generically develops the vacuum expectation value 



with Tr0 = 0. The gauge invariant order parameters Uk = (O|Tr0'^'|O) parameterize the Coulomb 
branch of the vacuum manifold. They define a scale in the theory with respect to which one 
can discuss the issue of a low energy effective action. This action takes into account one-loop 
perturbative correction as well as the whole instanton series and it is governed by the Riemann 
surface of the genus A"c — 1 for SU{Nc) gauge group. The same Riemann surface appears as the 
spectral curve of a classical integrable many-body system (see for a review) which turns out 
to be the Calogero type systems or some spin chains. The period matrix Tij[uk) of this Riemann 
surface yields the effective coupling constants of the gauge theory. 
For example, in the SU{2) case one has 




(51) 



4> = diag((/)i, ....,4>nJ , 



(52) 



+ 



e{u2) 



(53) 
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where g^{u2) and 9{u2) are the effective couphng and 6'— term respectively. The spectrum of the 
BPS states in this theory is given by 



Mnm = \na{u2) + maD{u2)\ (54) 
where a and an are the periods of some meromorphic differential Ag on the spectral curve S 



SW' r 'sw ' 

A Jb 



which in this case is a torus 

^sui2): y^ = {x'-K%x-U2). (56) 
The prepotential JF defining the low energy effective action is determined by the relations 

re.iu2) = ^ = ^, a^ = _ (57) 

with Tes{u2) given by (||). 

The connection with integrable system emerges when one studies the dependence of the 
prepotential on the fundamental scale A 

^^Pu,^0H, (58) 

where (3 is the one-loop beta-function of the gauge theory. This equation has another interpre- 
tation as evolution equation of an integrable system where H coincides with the Hamiltonian of 
this system and In A is the evolution time variable. The integrable system provides the natural 
explanation for the appearance of the meromorphic differential Asw, which turns out to be the 
action differential in the separated variables Ag^ = pdx. Let us emphasize that for SUSY YM 
case the classical integrable system is relevant and the meaning of the quantum system and the 
corresponding spectrum for SUSY YM remains an open question. It should involve the quan- 
tization of the vacuum moduli space since the hamiltonian in the dynamical system is nothing 
but H = U2 = (Tr$^). Simultaneously, this parameter serves as the coordinate on the moduli 
space of the complex structures of the Riemann surfaces whiich means that the quantization of 
the integrable system is related with the quantization of the effective d=2 gravity. 

The Riemann surface S degenerates at some points on the Coulumb branch on the moduli 
space. After the soft breaking of J\f=2 SUSY down to J\f=l SUSY these points correspond to the 
vacuum states of the Af=l gauge theory. Some massless states condense at these points leading 
to the formation of a mass gap and to confinement. In the SU{2) case, Eq. (^), the Af=l 
vacua correspond to the points U2 = ±A^ at which the monopoles or dyons become massless and 
condense. 

Let us now consider a particular theory, namely the superconformal Af=2 SUSY YM with 
Nf = 2Nc massless fundamental hypermultiplets 0. The corresponding integrable 

system is described by the spectral curve Sjv^ 

y' = PUx)-Ax'''^il-p'iT,i)), (59) 
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where p^ird) is some function of a coupling constant of the theory (see Eq. (^) below) and the 
polynomial Pn^ depends on the coordinates on the moduli space u= {u2, ...,U]\fJ 



iVc 



PnM = Yl l^i^) ^"^^"^ = 2^"^^ + ^2 + - + g^v. , (60) 

h=Q 

where go = 2, gi = and other are some known functions of u. Their explicit form is not 
important for our purposes. The Seiberg-Witten meromorphic differential on the curve (|59D is 
given by 

X^^ = pdx = hi{ijj /x^") dx ^ (61) 

where y = u — x'^^" /uo. 

From the point of view of integrable models the spectral curve (^) corresponds to a classical 
XXX Heisenberg spin chain of length with the spin zero at all sites (due to qi = 0) and 
parameter p related to the external magnetic field, or equivalently, to the twisted boundary 



conditions |[T0|1. 



Let us compare the spectral curve (|59D for the superconformal N'=2 SUSY YM with Nf = 2Nc 
with the spectral curve ( p3D for iV— reggeon compound states in multi-colour QCD. It is amusing 
to observe that these two curves coincide if we make the following identification 

• The number of the reggeons N = Nc] 

• The integrals of motion of multi-reggeon system are identified as the above mentioned 
functions qkiu) on the moduli space of the superconformal theory; 

• The coupling constant of the gauge theory should be such that p{tci) = 0. 

Under these three conditions the both theories fall into the same universality class. 

Let us clarify the meaning of the last condition. Since the superconformal theory has a 
vanishing /3— function one can assign a definite value to the bare coupling constant Tc\. The 
function p(r) entering ( |5UD at r = r^i is given by some combination of the modular forms |TP| , 
which can be determined from the duality properties of the superconformal SU{Nc) gauge theory. 
It turns out that the duality groups are different for odd and even Nc- The group for even Nc is 
generated by 

T: r^r+1, S : r^-l/r. (62) 
The group for odd Nc has the same T— transformation and 

S : T ^ —- ; ; — -— (63) 

(2cos(7r/2iVe))2r ^ ' 



S' = 1 {ST-y^ = 1 (64) 



subject to the constaints 

The explicit formulae for p(r) can be expressed in terms of weight four automorphic forms [|l 
for the duality group 

G,{t) + H,{t) 
P^"^ = GAr)-HAr) ' 
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where G4(r) is the Poincare series 

G,{r) = V(cr + d)-' , H,{r) = -^^G, ( -— ) , (66) 



c,d 

with a = 1 for even A''^ and a = 4 cos^(7r/(2A^c)) for odd A'^c- Here the summation goes over all 
pairs of integers (c, d) which occur in modular transformations r — (ar + 6)/ (cr + d) generated 
by T and STS~^. Finally, using ( |65D one finds that the condition p{tci) = is satisfied for 

1 Z TT 

Tel = — I — tan——. (67) 

The moduli of this superconformal theory are defined by the vector of the scalar condensates 
u and bare coupling constant Td. The spectral curve (|59D degenerates when p cxd or when the 
discriminant of the curve is zero. Let us consider simplest nontrivial case Nc = 3 which on the 
QCD side corresponds to the Odderon (A^ = 3) state. In this case the discriminant is given by 



2q? 
27q: 



Qs' = (6^ 



where / = 1 — p^. The curve degenerates at four values of among which three p^ = oo, and 1 
are universal and the last one depends on the moduli q2 and qs on the Coulomb branch. The first 
one is a strong coupling point corresponding to Td = 0. The second one is the so-called strong 
coupling orbifold point in which Xd is given by (|67D . It is this point on the moduli space which 
corresponds to the Regge limit of multi-colour QCD. Finally, the third one is a weak coupling 
point Tci — *■ oo. For A^^c = 3 the explicit expression for the function p('r), Eq. (|65D, is given in 
terms of the Dedekind function by fll 



"W-T^- /^M^^jn^wj '''' 

The strong coupling orbifold point p(rci) = describing the Odderon state in QCD occurs at 

Let us note this point becomes the Argyres-Douglas point when the moduli dependent singularity 
of the spectral curve p^(r) = 1— [1 -|- 2gf/(27g|)]^ collides with the strong coupling orbifold point 
([701). This happens when either 52/^3 = 0, or 



Going over to the general case of the SU{Nc) gauge theory one finds fl^, |l^ that indepen- 
dently on u, the spectral curve (|^) has three universal singularities at p^(rci) = cxd, and 1 
which have the same interpretation as in the case A"c = 3. 

Finally we would like to note that there is the following intriguing fact. In the case of A"c = 2 
which corresponds on the QCD side to the BFKL Pomeron state, the effective coupling constant 
is given in the weak coupling regime by the expression 

res = rd + ^^ + Vc,-e2^-'=-, (72) 
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where the second term is due to a finite one-loop correction I^TJ and the rest is the sum of 
instanton contributions. It is amusing that this one-loop correction to the coupling constant 



91e 



9l 



l + ^41n2 



+ 



(73) 



coincides with the expression for the intercept of the BFKL Pomeron P9| ) after one identifies 
bare coupling constant in the superconformal theory with the t'Hooft coupling constant in QCD. 

It would be very interesting to understand a proper QCD interpretation of the nonzero values 
of p in (|59|). From the point of view of spin chains it corresponds to the external magnetic field, 
or equivalently to the twisted boundary conditions |]TD|. The natural conjecture for p 7^ to 
occur is to consider the compound states in multi-colour QCD in which one of the reggeized 
gluons is replaced by the pair of reggeized quark and antiquark Similar picture takes place 



in the evolution equations for three-quark baryonic light-cone operators which will be discussed 
in Section 6. In that case we shall deal with a real SL{2, M) spin chain whose spectrum contains 
special eigenstates (with g = in Eq. (|109| ) below) for which a pair of quarks effectively behaves 
as one spin site @] . For generic values of p we can have a genus one curve even for a two particle 
system and the whole power of the duality group can be used. For example, the behaviour of the 
system near the strong coupling orbifold point p = can be related to the dual system near the 
weak coupling point p = 1- Certainly this possibility of having p 7^ in QCD deserves further 
investigation. 



4.2 Brane picture for the Regge limit 

Let us turn now to a stringy/brane picture for the Regge limit in multi-colour QCD. 

To warm up we would like to recall the brane description of the low-energy dynamics of the 
M = 2 SUSY YM. In the IIA framework the pure gauge theory is defined on the worldvolume of 
Nc D4 branes with the coordinates (xq, xi,X2, ^3, xq) stretched between two NS5 branes with the 
coordinates (xq, xi, 0:2, X3, 0:4, 2:5) and displaced along the coordinate xq by an amount inversely 
proportional to the coupling constant, 5xq = 1/0^ |2^. The coordinates at which the D4 branes 
intersect with the (X4, X5) complex plane define the vacuum expectation value of the scalar fields. 
This picture agrees perfectly with the RG behaviour of the coupling constant and yields the 
correct beta function in the gauge theory. The Riemann surface S discussed above describes the 
vacuum state of the theory and the spectrum of the stable BPS states. The lifting to the M 
theory picture leads to emergence of a single M5 brane with the worldvolume i?^ x S. 

In our case, we also need to incorporate into this picture branes corresponding to the fun- 
damental matter with Nf = 2Nc. There are two ways to do this: either using semi-infinite D4 
branes lifted into M5 brane in the M theory, or using D6 branes which induce the nontrivial KK 



monopole background for the M2 brane wrapped on the Riemann surface p3|] . As was shown 
in |jlO| in the latter case the resulting brane picture remains consistent with the integrable spin 
chain dynamics and we shall stick to this case. 

The explicit metric of the KK background in the M theory involving (X4, X5, Xg, Xio) coordi- 
nates has the multi-Taub-NUT form 

ds"^ = }^df^ + ^{dT + Adf)^ (74) 
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where r = (0:4, X5, Xe), r = xio and A is the Dirac monopole potential. The magnetic charge 
comes from the nontrivial twisting of bundle over R^. The function V behaves as 



V=l+y — — (75) 

i=l 

where = {x\,xI,Xq) are the positions of six-branes. For a superconformal case one must have 
5 = and positions of sixbranes in xq direction are irrelevant. 
Finally let us remind the brane interpretation of the BPS spectrum in the J\f = 2 SUSY YM. 
One way to realize these states is to consider the self-dual noncritical strings (which come from 
M2 branes ending on the surface) wrapped around the Riemann surface Since one has to 
keep some amount of SUSY the wrapping of the string should be along geodesies in some metric. 
This is equivalent to looking for the geodesies on u plane (in JIB approach), where u is the same 
as in Eg . (|30|) and is related to y = exp(— xe + ix^) and x = 0:4 + ix^ as y = u — x^^" juj. The 
metric on uj plane is defined in terms of Seiberg-Witten differential 

where Ag^ = pdx = (Ag^)*. Closed geodesies in the w— plane can be immediately lifted to 
closed curves on the Riemann surface. Hence the existence of the BPS state is the related to the 
existence of the closed geodesies in the (n, m)— cohomology class. Another important realization 
of BPS states arises from IIB/F theory picture. Namely the (n, m) BPS states come from the 
(n, m) strings which are stretched between the origin (which is generically split nonperturbatively 
into two strong coupling singularities) and the position of D3 brane we are living on [Q. Now 



the nontrivial metric on the Coulomb branch of the moduli space is involved. It has the form of 
the cosmic string metric with the conical singularity which can be brought into the flat form on 
the cover. In this picture SL{2, Z) invariance of the spectrum follows from the invariance of the 
IIB string theory. 

Let us turn now to our proposal for the brane realization of the Regge limit. We shall explore 
the brane representations for the Nf = 2Nc theory known in the IIA/M theory [^] and IIB/F 
theory p2[ . However unlike the SUSY case where the spectral curve is embedded in the internal 
"momentum" space the spectral curve of the noncompact spin chain is placed in the phase space 
involving both impact parameter plane as well as momenta. 

We shall first consider IIA/M type picture which is reminiscent to the realization of SYM 
theory via two NS5 and Nc D4 branes. We suggest that the coordinates involved in "HA" picture 
are the transverse impact parameter coordinates xi,X2 and rapidity A = In k^/k^. Transverse 
coordinates are analogue of (x4,X5) coordinates in SUSY case while rapidity substitutes the xq 
coordinate. 

Now let us make the next step and suggest that similar to SUSY case the single brane is 
wrapped around the spectral curve of XXX magnet and two "hadronic planes" together with 
N "Reggeonic strings" are just different projections of the single M2 brane with worldvolume 
i? X S. The coordinates involved into configurations are xi + ix2 and y = e~('*'+"'io) where xio 
is the "M-theory" coordinate. More precisely y can be identified with the Bloch-Floquet factor 
u in the definition of the Baker- Akhiezer function (|28|). The correspondence between "HA" and 
"M-theory" picture is shown in Fig. 3. 

Let us emphasize once again that the brane configuration for Regge limit contrary to SYM 
case partially involves the coordinate space. More precisely the geometry of Nf = 2Nc theory is 
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Figure 3: Spectral curve 



determined by the following parameter |23 



e = - ■ (77) 

Here A+ and A_ are asymptotic positions of five-branes defined by the large x behaviour of the 
curve 

iu oc A±x^^ . (78) 

A± can be found as roots of the equation 

Xl + X^ + \il-p') = (79) 

Since Regge limit corresponds to the strong coupling orbifold point, p = 0, the value of ^ is fixed 
as ^ = oo. This corresponds to the coinciding branes at infinity. 

Finally, the M theory brane picture for the Regge limit involves M5 brane corresponding to 
the vacuum state of the QCD. We can not say how it is placed precisely as the minimal surface 
in the internal seven dimensional space since the corresponding geometry is unknown yet. The 
new ingredient - M2 brane share the time direction with the M5 brane and wrapping around 
the Riemann surface which is embedded into two-dimensional complex "phase space" with the 
multi-Taub-NUT metric determined by KK monopoles with the magnetic charge 2N, which is a 
double number of reggeized gluons participating in the scattering process. 



5 Quantum spectrum and S-duality 

The 5— duality is a powerful symmetry in the SUSY YM theory which allows us to connect 
the weak and strong coupling regimes. The effective coupling in this theory coincides with the 
modular parameter of the spectral curve of the underlying classical integrable model. As a 
consequence, the 5— duality transformations in the gauge theory are translated into the modular 
transformations of the spectral curve describing complexified integrable system. A formulation 
of the 5— duality in the latter system naturally leads to an introduction of the notion of the dual 
action [|^. 

So far the S*— duality was well understood only for classical integrable models. In the case 
of multi-colour QCD in the Regge limit the situation is more complicated since the duality has 
to be formulated for a quantum integrable model. The integrals of motion take quantized set 
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of values and the coordinates on the moduh space are not continuous anymore. Therefore the 
question to be answered is whether it is possible to formulate some duality transformations at 
the quantum level. 

To study this question let us propose the WKB quantization conditions which are consistent 
with the duality properties of the complexified dynamical system whose solution to the classical 
equations of motion are described by the Riemann surface. We recall that the standard WKB 
quantization conditions involve the real slices of the spectral curve 

pdx = 2nh{ni + 1/2) (80) 

where rii are integers and the cycles Ai correspond to classically allowed trajectories on the phase 
space of the system. In our case the coordinate x is complex and arbitary point on the Riemann 
surface is classically allowed. As a result the general classical motion involves both A— and 
5— cycles on the Riemann surface. This leads to the following generalized WKB quantization 
conditions (see also pO[| ) 



Re ® pdx = nhrii , Ke (b pdx = nhmi , (81) 



where the "action" differential was defined in (pll) . Note that in the context of the SUSY YM 
this condition would correspond to the nontrivial constraints on the periods (|55| ) and on the mass 
spectrum of the BPS particles (|^). It is clear that the WKB conditions (|8l| ) imply the duality 

Ai ^ Bi and rii ^ uii. 

The conditions (^Tj) are analogous to the WKB quantization for a particle in a multi-dimensional 
case. We have exactly this situation since we are dealing with complexified phase space. For a 
particle in higher dimensions one considers the representation of the wave function in terms of 
the multivalued functions Ak and Sk on the x-space 

m{x) = Y, M^) exp {\Sk{x)\ . (82) 

k ^ ' 

Here Ak{x) and Sk{x) are different branches of the multivalued functions A{x) and S{x). To 
obtain the WKB quantization conditions one considers the covering space on which A{x) and 
S{x) are single valued and requires that the wave function \l/(x) should resume its original value 
after encircling along all noncontructable cycles Ci on the covering manifold. This leads to 

p-dx+(p d\n A{x) = 2nhni . (83) 

The second integral is equal to the sum of two terms related to the topological invariants of 
the classical trajectory: the number of intersections of Ci with hypersurface of the caustics 
and the number of intersections of Ci with the hypersurface of the turning points. In our case 
p ■ dx = 2 Re p dx and §^ d In A{x) is an even integer. 



Let us consider the quantization conditions (q1|) in the simple case of the Odderon N = 3 
system. The spectral curve (p3| ) is a torus 

= {2x^ + q2X + qsY - 4a;^ (84) 
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where q2 is given by (|35D and (^6]) while is the complex integral of motion to be quantized. 
The quantization conditions (|8T| ) read (for h = 1) 



Re 0(53) = vrn, Rea£)(g3) = vrm (85) 

where n and m are integer. These equations can be solved for large values of g'l/q'l ^ 1, for 
which the expressions for the periods 0(^3) and 00(^3) are simplified considerably. The explicit 
evaluation of the integrals ( p5D in this limit yields 



Substituting these expressions into (|^) one finds 




1/3 _ r3(2/3) 



= 7^ + • (87) 



where £1 = n and £2 = ?^ — 2m. The same quantization conditions can be also written as 

(4 + 3^2) ■ a{q-i) - U2 ■ aDiqs) = vr(^? + 3il) . (88) 

The WKB expressions ( |87| ) are in a good agreement with the exact expressions for quantized ^3 
obtained from the solutions of the Baxter equations in |20[ . 



There is a simple relation between the modular parameter of the curve with the periods a 
and ao 

an = TefT a . (89) 

It follows from the fact that a and are the solutions to the Picard-Fuchs equations which in 
the Odderon case is the second order differential equation. One can show that the Wronsian for 
this equation aa'^ — a^a' vanishes and using t^s = a'^/a' one gets (|89D. Due to this fact half of 
quantization conditions can be formulated in terms Tes instead of qd- 

Let us emphasize that the point at the moduli space corresponding to the degeneration of the 
torus for the Odderon case does not appear in the quantum spectrum. From the point of view 
of the gauge theory this means that the appearance of the massless states is forbidden. 

In the general multi-reggeon case we have to consider the quantization conditions (|8TD on 
the Riemann surface (|3^) of the genus (A^ — 2) which has the same number of the A— and 
S— cycles. In result the spectrum of the integrals of motion q^, q^ is parametrized by two 
[N — 2)— component vectors n and m. In the SUSY YM case these vectors define the electric 
and magnetic charges of the BPS states. In the Regge case the physical interpretation of n and 
rh is much less evident. Let us first compare the electric quantum numbers in the two cases. In 
Regge case it corresponds to rotation in the coordinate space around the ends of the Reggeons. 
This picture fits perfectly with the interpretation of the electric charge in SUSY YM case. Indeed 
VEVs of the complex scalar take values on the complex plane which is the counterpart of the 
impact parameter plane and the rotation of the phase of the complex scalar is indeed the "electric 
rotation" . 

To get some guess concerning the "magnetic" quantum numbers it is instructive to check 
the geometrical picture behind them in the simplified "IIA" picture. All states corresponding 
to the "electric" degrees of freedom are related to fundamental strings encircling "reggeoinic" 
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tubes (see Fig. 3) and don't feel the hadronic planes. However the "magnetic" states as is well 
known from SUSY YM case are represented by the membrane stretched between two strings and 
two hadronic planes. Therefore these states are sensitive to hadronic quantum numbers. More 
detailed interpretation of "magnetic" degrees of freedom in the Regge regime has to be recovered. 

Let us emphasize that due to the stringy realization of the abovementioned BPS spectrum 
we obtain rather transparent picture of the WKB quantization in the stringy terms. Remind 
that, in SUSY YM the mass of BPS state Mn^m is given by the energy of the dyonic (n, m) string 
stretched between the value of U2 corresponding to a given VEV of a scalar field and the point 
^2 = where orientifold and all D7 branes are placed. Hence, the WKB quantization at least in 
the Odderon case can be formulated as the quantization of the energy of the dyonic strings at 
the complex gs plane and the spectrum inherited the duality from S-duality of HB string. 



6 Stringy /brane picture and the calculation of the anoma- 
lous dimensions 

In the previous Sections we have demonstrated that integrability properties of the Schrodinger 
equation for the compound state of Reggeized gluons give rise to the stringy/brane picture for 
the Regge limit in multi-colour QCD. There is another limit in which QCD exhibit remarkable 
properties of integrability. It has to do with the scaling dependence of the structure functions of 
deep inelastic scattering and hadronic light-cone wave functions in QCD. In the both cases, the 
problem can be studied using the Operator Product Expansion and it can be reformulated as a 
problem of calculating the anomalous dimensions of the composite operators of a definite twist. 
The operators of the lowest twist have the following general form 

0%{^) = (i//^)'=^$i(0)(i/Z})'=^$2(0)(yD)^-'=-'=^$3(0), (90) 

where k = {ki, ^2) denotes the set of integer indices ki, is a light-cone vector such that = 0. 
$fc denotes elementary fields in the underlying gauge theory and = — iA^ is a covariant 
derivative. The operators of a definite twist mix under renormalization with each other. In 
order to find their scaling dependence one has to diagonalize the corresponding matrix of the 
anomalous dimension and construct linear combination of such operators, the so-called conformal 



operators p4| , [45 



o^^::;{o) = j2Ck,, -0^,^(0). (91) 

k 

A unique feature of these operators is that they have an autonomous RG evolution 

Here is a UV cut-off and 'y^^g is the corresponding anomalous dimension depending on some 
set of quantum numbers q to be specified below. It turns out that the problem of calculating the 
spectrum of the anomalous dimensions '-fN,q to one-loop accuracy becomes equivalent to solving 
the Schrodinger equation for the S'L(2,M) Heisenberg spin magnet. The number of sites in the 
magnet is equal to the number of fields entering into the operators under consideration. 
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To explain this correspondence it becomes convenient to introduce nonlocal light-cone oper- 
ators 

F{zi, Z2) = <^i{ziy)<^2iz2y) , F{zi, Z2, Z3) = ^i{ziy)<^2iz2y)^3{z3y) . (93) 

Here is a light-like vector (y^ = 0) defining certain direction on the light-cone and the scalar 
variables Zi serve as a coordinates of the fields along this direction. The fields ^i{ziy) are 
transformed under the gauge transformations and it is tacitly assumed that the gauge invariance 
of the nonlocal operators F{zi) is restored by including the Wilson lines between the fields in 
the appropriate (fundamental or adjoint) representation. The conformal operators appear in the 
OPE expansion of the nonlocal operators ( P^D for small zi — Z3 and 22 — ^3- 

The field operators entering the definition of F{zi) are located on the light-cone. This leads 
to the appearance of the additional light-cone singularities. They modify the renormalization 
properties of the nonlocal light-cone operators (^) and lead to nontrivial evolution equations 
which as we will show below become related to integrable chain models. We notice that there 
exists the following relation between the conformal three-particle operators (|9TD and the nonlocal 
operators (0) 

0^-/(0) = vl/^,,(a,„ a,3)F(zi, Z2, zs) , (94) 

Zi=0 

where \l/7v,g(a;i, 0:2, X3) is a homogenous polynomial in Xi of degree N 

^7v,,(xi, X2, xs) = Cfc,<, • x'l'x'2'x^-''-''' (95) 
fc 

with the expansion coefficients Ck,q defined in (pT]). Similar relations hold for the twist-2 op- 
erators. The problem of defining the conformal operators is reduced to finding the polynomial 
coefficient functions '^N,q{xi) and the corresponding anomalous dimensions 7Ar,ij. Using the renor- 
malization properties of the nonlocal light-cone operators (|93|) one can show that to the one- 
loop accuracy the QCD evolution equation for the conformal operators (^) can be rewritten in 
the form of a Schrodinger equation 

n ■ ^N,gix,i) = 7iV,g^JV,g(xO , (96) 

where the Hamiltonian Ti. acts on the Xj— variables which are conjugated to the derivatives dz^ 
and, therefore, have the meaning of light-cone projection {y ■ pi) of the momenta pi carried by 
particles described by fields ^{ziy). 

The Hamiltonian Ti has the following remarkable properties. 

(i) This Hamiltonian is a sum of two-particle Hamiltonians 

where Hij acts on the coordinates of particles i and j. 

(ii) All two-particle Hamiltonians Hij are invariant under SL{2, R) transformations 

[Hij,Si + Sj]=0 (98) 
where one-particle SL{2,R) generators Si act on the space of functions \l/(xj) as 

Skfi'^ixk) = {xkdk+ ik)'^{xk) , 
S'fc,+^(xfc) = -Xk'^{xk) , 

Sk,--^{xk) = [xkdl + 2jkdk)-^{xk). (99) 
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Let us note that these generators where Xk play the role of momenta variables are dual to gen- 
erators (lU) where Zk were coordinates. This symmetry of one-loop evolution kernels Hij follows 
from the invariance of classical QCD Lagrangian under the group of conformal transformations 
which is reduced to its SL{2, R) subgroup on the light-cone 



az + b 



cz + d 

with ad — be = 1 and a, b, c, d real. Here 



(100) 



jk = ^{dk + Sk), (101) 

where dk is a canonical dimension of a field ^k{zy) and Sk is a light-cone component of it spin. 
For example, dq = 3/2 and Sq = 1/2 and jq = 1 for quarks. Due to (^) Hamiltonians Hij are 
functions of the two-particle Casimir operators 

{S, + S2y = Ju{Ji2-l). (102) 
For example when $i and $2 are quark fields of the same chirality 

F^p{zuZ2) = J2{qj^Uzw)Uql)p{z2y) (103) 

i=l 

with qj = (1 -|- 75)^4/2, the two-particle Hamiltonian is given by p5| 

Hl2 = -Cf [HqqiJu) + 1/4] , Hqq{J,2) = ^(Ju) " V^(2). (104) 

where Cp = {N^ — l)/{2Nc). The eigenfunctions for this Hamiltonian are the highest weights of 
the discrete series representation of the SL{2, R) group 

¥^\x,,X2) = ix, + X2fcT f^l^) (105) 

\Xi + X2 / 

where C^"^ are Gegenbauer polynomials. The corresponding eigenvalues define the anomalous 
dimensions of the twist-2 mesonic operators built from two quarks with the same helicity 

A' 



7?^ = -CpmN + 2) - V^(2) + 1/4] = ^Cf 



^ 1 1 

.k=l 



(106) 



At large this expression has well-known asymptotic behaviour ~ OgCiT'/vrln A^. 

It is conformal symmetry which dictates that the two-particle Hamiltonian is a function 
of the Casimir operator of the 5'L(2,R) group, but it does not fix this function. The fact 
that this function turns out to be the Euler T/^-function leads to a hidden integrability of the 
evolution equations for anomalous dimensions of baryonic operators in very much the same way 
as emergence of the Euler ?/^-functions in the BFKL kernel (0) leads to the integrability of multi- 
colour QCD in the Regge limit. Namely, for baryonic operator built from three quark fields of 
the same chirality 

Fap^{zuZ2,Z3) = ^ eijkUqj)a{ziy){^q])fs{z2y){^ql)^{z3y) (107) 

i,j,k=l 
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the evolution kernel is given by @, ^ 



^(3) = ^ {(1 + i/N,) [HggiJu) + H,,{J23) + H,,{J3,)] + 3Cf/2} (108) 
Zn 



with Hqq given by ( |104| ). The Schrodinger equation ( p6D with the Hamiltonian defined in this 



way has a hidden integral of motion 

q = i {d^^ - d^^) (9^2 - 9^3) ((9^3 - X1X2X3 (109) 

and, therefore, it is completely integrable. This operator is dual to the operator gs for the 
N = 3— reggeon states, Eq. ([T3|), in the same fashion as the SL{2) generators were related in the 
two cases. Similar to the Regge case, one can identify (|10(j|) as the Hamiltonian of a quantum 



XXX Heisenberg magnet of SL{2, M) spin jq = 1. The number of sites is equal to the number of 
quarks. 

Based on this identification we shall argue now that the calculation of the anomalous di- 
mensions can be formulated entirely in terms of Riemann surfaces which in turn leads to a 
stringy/brane picture. It is important to stress here the key difference between Regge and light- 
cone limits of QCD. In the first case the impact parameter space provides the complex plane 
for the Reggeon coordinates and we are dealing with a (2 + 1)— dimensional dynamical system. 
In the second case the QCD evolution occurs along the light-cone direction and is described by 
a (1 + 1)— dimensional dynamical system. As a consequence, in these two cases we have two 
different integrable magnets: the 5'L(2,C) magnet for the Regge limit and the SL{2,'R) magnet 
for the light-cone limit. The evolution parameters ("time" in the dynamical models) are also dif- 
ferent: the rapidity In s for the Regge case and the RG scale In fi for the anomalous dimensions 
of the conformal operators. 

Our approach to calculation of the anomalous dimensions via Riemann surfaces looks as 
follows. For concreteness, we shall concentrate on the evolution kernel (|108|) . Similarly to the 
Regge case, one starts with the finite-gap solution to the classical equation of motion of the 
underlying SL{2,M.) spin chain and identifies the corresponding Riemann surface 

uj- — =2x^-(N + 2)(N + 3)x + q, u = x^ (110) 

UJ 

where q is the above mentioned integral of motion ( |109| ) and N is the total SL{2,W) spin of 
the magnet, or equivalently the number of derivatives entering the definition of the conformal 
operator (|9T|) . In distinction with the SL{2, C) case, Eq. (^61), the total spin takes nonnegative 
integer values. Note that the Riemann surface corresponding to the three-quark operator has 
genus g = 1, while g = for the twist 2 operators. 

At the next step we quantize the Riemann surface as follows. We replace p = id/dx and 
impose the equation of the complex curve as the operator annihilating the Baxter function 

(e^9/9^ + e-*^/^^')x2Q(x) = [2x^ -{N + 2){N + 3)x + q]Q{x). (Ill) 

This leads to the Baxter equation for the Heisenberg SL{2, M) magnet of spin j = 1 

{x + ifQ{x + i) + {x- ifQix -i) = [2x^ -{N + 2){N + 3)x + q] Q{x) . (112) 

Similar to the Baxter equation in the Regge case, this equation does not have a unique solu- 
tion. To avoid this ambiguity one has to impose the additional conditions that Q{x) should 
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be polynomial in x. This requirement leads to the quantization of the integral of motion. The 
resulting polynomial solution Q = Qq{x) has the meaning of the one-particle wave function in 
the separated variables x which in the case of the SL{2,M.) magnet take arbitrary real values. 



Given the polynomial solution to the Baxter equation (|112|) , one can determine the eigen- 
spectrum of the Hamiltonian ( p,08| ) and in result the anomalous dimensions of the corresponding 
baryon operators 



(3) 

rN,q 
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[(1 + l/N,)£M,g + 3Cf/2] , £ 



N,q 



— I- 



.Q'a' 



(113) 



parameterized by the eigenvalues of the integral of motion ( |109| ) given by 

■Qq{i) - Qq{-i) 



q = —I- 



QM 



;ii4) 



The solution of the Baxter equation simplifies greatly in the quasiclassical approximation 
which is controlled by the total S'L(2,R) spin of the system A^. For N ^ 1 the spectrum of the 

im. 



integral of motion q is determined by the WKB quantization condition 



pdx = 2'K{n + 1/2) + 0{l/N) 



(115) 



where p was introduced in ( p. 1 1| ) . Here integration goes over the A-cycle on the Riemann surface 
defined by the spectral curve ( p.ll|) . This cycle encircles the interval on the real x— axis on which 



e^' I > 1. Solving ( |115| ) one gets 



ATS 



<1 



(116) 



Comparing this expression with (fflD and taking into account that in this case q2 = — (A^ + 3)(A^ + 
2) and gs = g we conclude that for N (yo the system is approaching the Argyres-Douglas point. 
Note also that the WKB quantization conditions ( |115| ) involve only the A— cycle on the Riemann 
surface and unlike the Regge case there is no S"— duality in the quantum spectrum. 

Finally, the spectrum of the anomalous dimensions in the WKB approximation is given by [|Tl|, 
Si 

3 

SN,q = 2 In 2 - 6 + 67e + 2Re ^ ^/^(l + i6k) + 0{N~^), (117) 

fc=i 

where 5k are defined as roots of the following cubic equation: 



26l-{N + 2){N + 3)6k + q 







:ii8) 



and q satisfies ( [1161 ). 

What can we learn about stringy picture from this information about anomalous dimensions? 
Let us remind that in the spirit of string/gauge fields correspondence |^ the anomalous di- 
mensions of gauge field theory operators coincide with excitation energies of a string in some 
particular background. An important lesson that we have learned from the analysis of the the 
two- and three-quark operators is that in the first case the anomalous dimensions are uniquely 
specified by a single parameter A^ which define the total SL{2,'R) spin. In the second case, a 
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new quantum number emerges due to the fact that the corresponding dynamical system is com- 
pletely integrable. The additional symmetry can be attributed to the operator q defined in (|109|) . 
From the point of view of a classical dynamics this operator generates the winding of a particle 
around the A— cycles on the spectral curve. Within the string/gauge fields correspondence one 
expects to reproduce these properties using a description in terms of the same string propagating 
in different backgrounds. One is tempting to suggest that different properties of the anomalous 
dimensions of the two- and three-particle opereators should be attributed to different properties 
of the background. In the case of the twist-2 the anomalous dimensions depend on integer 
which in the classical system has an interpretation of the total S'L(2,M) angular momentum of 
the system. On the stringy side the same parameter has a natural interpretation as a string 
angular momentum. 

When this paper was in a process of preparation the interesting paper concerning the stringy 
derivation of the anomalous dimensions of the twist-2 operators [0 appeared.^ It was shown 



that the classical solution to the equation of motion in the cx-model on AdS^^ x 5*^ provides the 
anomalous dimension of the twist-2 operators with the large spin A^. Since our derivation of the 
same anomalous dimension is seemingly different it is very instructive to compare two approaches. 



The logic behind derivation in implies that one considers the stringy a model perturbed by 
some vertex operator representing the operator under consideration on the gauge theory side. 
Then the energy on the solution to the classical equations of motion in a— model involving the 
radial coordinate in the bulk provides the anomalous dimension of the operator. 
In our approach we have the following correspondence 

operator Riemann surface 

twist of the operator ■^=r- genus of the Riemann surface 
calculation of the anomalous dimension quantization of the Riemann surface 

It seems that the Riemann surfaces whose moduli (the integrals of motion of the spin chain) 
define the anomalous dimensions of the corresponding operators describe the a— model solutions 
found in H^. The precise relation between two approaches need to be clarified further and will 



be discussed elsewhere. 

As we have seen the quantization of the Riemann surface can be performed most effectively 
in terms of the Baxter equation. It is worth noting that the solution to the Baxter equation 
can be identified as a wave function of DO brane |4^. Quantization conditions arise from the 



requirement for the wave function the DO brane probe in the background of the Riemann surface 
to be well defined. 

In the case of the light-cone composite operators we have to incorporate into the stringy 
picture a new quantum number which is parameterized by an integer n, Eq. ( p.l6| ), i.e. string 
excitation spectrum has now different sectors parameterized by this integer. The natural way 
to interpret these sectors is to identify n with the winding number of a closed string. The 
corresponding background for such scenario is offered by the Riemann surface itself with the 
string wrapped around the A— cycle. It is an interesting open question if one can interpret a 
momentum in WKB quantization condition ( |115|) as a momentum of a string T-dual to the 
string with the winding number n. 



^Note that some attempts to develop the cr-model representation for high energy QCD were undertaken a long 
time ago 
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Since the spectrum of the anomalous dimensions in QCD coincides with the spectrum of the 
SL[2,M.) spin chain Hamiltonian it would be interesting to explore further the symbolic relation 

-^^string OC H^pi^ (119) 

where string propagates in the background determined by the Riemann surface of the spin chain. 
It is known that hamiltonian formulation of the spin chains is closely related to the Chern- 
Simons (CS) theory with the inserted Wilson lines. The number of the sites in the spin chain 
corresponds to the number of the Wilson lines. The gauge group in the CS theory is the 
symmetry group of the magnet. Here it is the SL{2, M) and group manifold is an AdS^ space. 
Hence the spin chain describes the motion of points in AdSs space. 



7 On the dual bulk representation of the reggeon 

Let us make some comments about a dual gravity description of multi-colour QCD in the Regge 
limit. We can not present the complete solution of the problem due to lack of the proper 
(super) gravity background for nonsupersymmetric theories but mention some features which 
seem to be important for the behaviour of the scattering amplitudes. First of all we want to 
emphasize that in spite of the large Nc limit one has to deal with Riemann surfaces whose genus is 
determined by the number of reggeons. In the scattering process the relevant geometry is captured 
by the M2 brane wrapped around the Riemann surface in the background Taub-NUT metric. 
Two transverse coordinates and two momenta are involved in this Taub-NUT background. It is 
this part of the total metric which provides the intercepts of the multi-reggeon states. We expect 
that back reaction of "scattering" M2 brane on the "vacuum" M5 brane can be neglected. 

The M2 brane modifies the background and it is known that its near horizon geometry involves 
AdS4^ factor. One of the M2 coordinates is the time direction. Because here we are dealing with 
static properties of the system one can study them at any moment in time. As a consequence 
we are interested in the fixed time near horizon geometry which gives us AdS^ space. This could 
help to determine a proper CFT on the transverse plane but an exact respective conformal 



field theory is still unknown and CFT interpretation of reggeized gluons and their compound 
states in pomeron and odderon sectors remains obscure. 

Let us formulate here a conjecture that in an AdS/CFT description reggeized gluons are 



described by AdS^ singletons [^]. This conjecture is based on the observation that the correlation 
function of two reggeized gluons are logarithmic in transverse space and they can be identified 
with zero dimension logarithmic operators. It was shown some time ago that in the AdS/CFT 
correspondence the logarithmic operators on the boundary correspond to the singletons in the 

bulk ig. 



To see the logarithmic correlator to appear let us consider the amplitude for the scattering 
of two quarks with infinite energy 

W^/ = {0\TWl'{0)W!^^{z)\0) . (120) 

Here, Wilson lines and are evaluated along infinite lines in the direction of the quark 
velocities fi and V2, respectively: 



W+{0) = Pexp (^i j davi- A{via)^ 
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W^z) = Pexp i^t j d(3v2- A{v2p + z)j , (121) 

and the integration paths are separated by the impact vector z = {0'^,0~,z) in the transverse 
direction. 

To understand the 2;— dependence, let us consider as an example the one-loop calculation of 
W in the Feynman gauge. One gets 

Wi.ioop = i(g)i + ie(g) e)^^T{D/2 - i)A^-^ 

^ da r ^-^1^, - (122) 

In this expression, gluon is attached to both Wilson lines at points via and V2l3 + z and we 
integrate the gluon propagator V1V2D ^y{via—V2[3—z) over positions of these points. To regularize 
IR divergences we introduced the dimensional regularization with Z) = 4 + 25, (e > 0) and A 
being the IR renormalization parameter. There is another contribution to Wi^ioop corresponding 
to the case when gluon is attached by both ends to one of the Wilson lines. A careful treatment 
shows that this contribution vanishes. Performing integration over parameters a and (3 we get 

W^i-;oop(7,A'z') =/»/+ (r ®r) — (-i7rcoth7)r(£)(7rA2i'2)-^ (123) 

TT 

where (^1^2) = cosh 7. 

The integral over a and f3 in ( [122| ) has an infrared divergence coming from large a and /?. In 
the dimensional regularization, this divergence appears in Wi-ioop as a pole in {D — 4) with the 
renormalization parameter A having a sense of an IR cutoff. We see that in an octet channel we 
got a logariphmic correlation which depends only on two transverse coordinates 

(r ®r) — (i7rcoth7)ln(7rA2i'2) (124) 

71 

which allows us to make a conjecture that reggeized gluon is described by some logarithmic 
operator with zero dimension which correspond to a singleton field in AdS^. 

The arguments above suggest that the CFT at the transverse plane could be logarithmic 
P^ . Also it is tempting to relate the SL(2,C) symmetry of the AdS^ space with the SL{2,C) 



symmetry of the original XXX spin chain but at this stage we do not have reasonable arguments 
in favor of this connection. Certainly these important issues deserve further investigation. 



8 Conclusions 

Let us summarize the obtained results. In this paper we proposed the stringy picture for multi- 
colour QCD in the Regge limit and on the light-cone in which hidden integrability plays the 
central role. Our approach is similar in many respects to the one which is known to be very 
successful in the description of the low energy limit of the A/'=2 SUSY YM theories. In the 
Regge limit we encounter a new situation in comparison with a SUSY case. One has to develop 
the quantum picture involving the dynamics on a Riemann surface as the quasiclassical limit. We 
argued that the whole configuration relevant in the Regge limit follows from M2 brane wrapped 
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around the spectral curve of the integrable system. Within our approach the Regge hmit turns 
out to be in the same universahty class as M=2 superconformal SUSY YM with Nj = 2Nc at 
the strong coupling regime. In both cases there is no natural place for the Aqcd type parameter: 
in SUSY case theory is conformal while in the Regge case we are dealing only with perturbative 
regime in generalized leading logarithmic approximation. 

The WKB quantization conditions providing the compound multi-reggion spectrum and the 
spectrum of anomalous dimensions for composite light-cone operators are formulated in the 
stringy terms. They imply that energy of strings properly wrapped around the spectral curve 
has to be quantized. Note that since the quantization of the integrable system amounts to the 
quantization of the moduli space of the complex structures of the Riemann surface quantum 
gravity in two dimensions is involved. The correspondence operator <;=^ Riemann surface plays 
an important role in our approach and is some new step toward a complete formulation of the 
gauge/string correspondence. 

One more issue which certainly deserves further investigation is the meaning of the unitary 
Froissart boundary. Let us emphasize that to get the scattering amplitude one has to sum over 
all Riemann surfaces of arbitrary genus, corresponding to the different number of Reggeons N 

A{s,t) = J2^N{s,t)= Yl M^^t) (125) 

N>2 gcnus=N-2 

in perfect agreement with the string picture. Hence, if the picture suggested for the Froissart 
saturation in is correct it would mean that the process of the black hole creation in the s- 
channel is equivalent to the closed string propagation in nontrivial background in the t-channel. 

Let us mention a few other important questions which remain open. One of them is to find 
more clear physical identification of the M-theory 10-th dimension involving the geometry of the 
spectral curve. At a moment we can accurately introduce it only via the Bloch-Floquet factor 
of the Baker- Akhiezer function but more physical interpretation is highly desirable. The related 
question concerns the interpretation of the "magnetic" quantum numbers amounted from the 
spectral curves. It would help in a wider application of the dualities in the Regge limit and 
possibly could shed the additional light on the s — t duality in the high energy amplitudes. 

We are grateful to S. Derkachov, Yu. Dokshitzer, A. Kaidalov, A. Leonidov, A. Manashov, 
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by grants CRDF-01-2108 and INTAS-00-334, I.K. is supported in part by PPARC rolling grant 
PPA/G/O/1998/00567 and EC TMR grant HPRN-CT-1999-00161, both I.K. and O.K. are par- 
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